REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and  maintaining  the 
data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing 
this  burden  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202- 
4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently 
valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE  3.  DATES  COVERED  (From  -  To) 

01-10-2006  Final  report  05-15-2005  -  11-14-2005 


4.  TITLE  AND  SUBTITLE  5a.  CONTRACT  NUMBER 

Hybrid  Molecomputer  Using  Vapor  Phase  Assembly  NA 

5b.  GRANT  NUMBER 

FA9550-05-1-0340 

5c.  PROGRAM  ELEMENT  NUMBER 

NA 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES)  8.  PERFORMING  ORGANIZATION  REPORT 

NUMBER 

Yale  University  NA 

Grant  and  Contract  Admin . 

155  Whitney  Ave  Suite  214 
New  Haven,  CT  06520-8337 


6.  AUTHOR(S)  5d.  PROJECT  NUMBER 

Reed,  Mark  A.  NA 

5e.  TASK  NUMBER 

NA 

5f.  WORK  UNIT  NUMBER 

NA 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

AFOSR 

4015  Wilson  Blvd 
Room  713 

Arlington,  VA  22203-1954 


12.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Public  distribution/availability 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 

AFOSR 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

AFRL- SR- AR- 06-0233 


14.  ABSTRACT 

Current  work  on  molecular  electronics  relies  almost  entirely  on  liquid  phase  processing.  Gas 
phase  processing  offers  far  superior  uniformity  and  reproducibility.  The  scope  of  this  work 
was  to  investigate  vapor  phase  process  that  will  fully  integrate  the  in-situ  steps  of  sample 
preparation,  active  organic  molecule  deposition,  and  gentle  gas -moderated  metallic  over¬ 
coating,  which  will  enable  our  proposed  multi-layer  stacked  device.  The  goal  of  this  work  is 
to  mate  these  active  molecular  gates  with  a  MOS-based  post -processing  circuit. 


15.  SUBJECT  TERMS 

molecular  electronics,  VPE 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

OF  ABSTRACT 

18.  NUMBER 
OF  PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

Dr.  Anne  Matsuura 

a.  REPORT 

none 

b.  ABSTRACT 

none 

c.  THIS  PAGE 

none 

NA 

19b.  TELEPHONE  NUMBER  (include  area 
code) 

703-696-6204 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39.18 


1 


Final  report  FA9550-05-1-0395 

Objectives 

The  objectives  of  this  effort  are  to  explore  the  use  of  a  functionalized  semiconductor  nanowire 
matrix  with  nucleic  acid  reporters  for  chem/bio  sensors.  The  milestones  were: 

•  demonstrate  unambiguous  detection  of  an  analyte  in  an  aqueous  environment  by  gating 
of  a  semiconductor  nanowire 

•  demonstrate  identification  of  an  analyte  using  a  functionalized  nanowire  process 

•  demonstrate  nucleotide  functionalization  of  a  semiconductor  nanowire. 

•  downselect  the  specific  nanowires  and  target  functionalizations 

Status 


All  milestones  listed  above  were  achieved.  The  program  is  continuing  into  Phase  II  funding. 

The  progress  is  detailed  in  the  next  section 

Accomplishments/New  Findings 

Background 

In  recent  years,  investigators  have  made  tremendous  strides  in  demonstrating  the  basic 
capabilities  of  molecular-based  electronics.  However,  this  work  consists  largely  of  simple 
circuits,  fabricated  under  carefully  controlled  laboratory  conditions,  by  extremely  well  trained 
investigators.  The  goal  of  this  program  is  to  develop  a  hybrid  molecular/solid-state  circuit 
configuration  that  will  be  fabricated  by  a  more  scaleable  and  production-worthy  vapor  phase 
fabrication  process.  The  hybrid  approach  builds  on  and  combines  existing  capabilities  of 
molecular  and  solid-state  circuits:  The  functional  unit  is  a  bistable  molecular  gate  consisting  of 
either  a  stacked  negative  differential  resistive  (NDR)  or  memory  molecular  device.  This  device 
has  already  been  successfully  demonstrated.  These  active  molecular  gates  will  then  be  mated 
with  a  MOS-based  post-processing  circuit  to  make  dense,  high-performance  cells.  This  union 
will  be  facilitated  by  new  vapor  phase  deposition  process  will  allow  vastly  more  controlled 
fabrication  of  the  fully  stacked  molecular  gates  upon  a  MOS-based  substrate  template. 

Current  work  on  molecular  electronics  relies,  almost  entirely,  on  liquid  phase  processing. 
Liquid  phase  processing  is  simple,  it  is  inexpensive,  and  it  builds  upon  a  large  body  of 
experience  with  self-assembled  molecular  monolayers  (SAMs).  We  now  want  to  drive 
molecular  circuits  to  a  complexity  that  will  match  (and  eventually  exceed)  that  of  solid-state 
microelectronics.  Early  microelectronic  processing  was  also  based  on  the  use  of  liquids. 
However,  as  circuit  complexity  and  area  increased,  liquid-based  processing  could  not  deliver  the 
uniformity  and  reproducibility  required  for  successful  fabrication.  Thus,  virtually  all  modem 
semiconductor  tools  are  based  on  the  use  of  gases.  Gas  phase  processing  offers  far  superior 
uniformity  and  reproducibility.  Indeed,  when  surface  reactions  dominate,  it  is  possible  to 
simultaneously  deposit  or  etch  the  surfaces  of  dozens  of  stacked  wafers.  More  subtly,  gases  (and 
ions)  allow  for  very  precise  and  shallow  etching.  This  is  essential  in  the  fabrication  of  a  circuit 
of  microprocessor-level  complexity:  microprocessors  cannot  be  fabricated  in  a  single  deposition 
step.  The  operation  of  the  device  cell,  and  interconnection  of  cells,  inevitably  requires  a  complex 
multi-level  structure.  It  is  thus  essential  that  the  technology  allow  for  repeated  reinsertion  into 
deposition  systems.  For  such  a  sequence  of  depositions  to  work,  one  must  be  able  to  clean  and 
prepare  underlying  structures,  in-situ,  for  the  deposition  of  new  materials.  Controlled  deposition 
and  etching  are  thus  essential. 

This  project  is  a  collaboration  between  University  of  Virginia  and  Yale  University.  Yale’s 
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responsibility  was  the  fabrication  and  characterization  of  molecular  devices;  UVA  had  the 
responsibility  of  vapor  phase  molecular  deposition.  This  report  focuses  on  the  fabrication  and 
transport  data. 

Progress  and  Accomplishments 

The  suggestion  and  demonstration”  of  utilizing  molecules  as  the  active  region  of  electronic 
devices  has  recently  generated  considerable  interest  in  both  the  basic  transport  physics  and 
potential  technological  applications  of  “molecular  electronics”.3,4  However  some  reports  of 
molecular  mechanisms  in  electronic  devices5'6  have  been  shown  to  be  premature  and  due  to 
filamentary  conduction,7  highlighting  the  fabrication  sensitivity  of  molecular  structures  and  the 
need  to  institute  reliable  controls  and  methods  to  validate  true  molecular  transport.  A  related 
problem  is  the  characterization  of  molecules  in  the  active  device  structure,  including  their 
configuration,  bonding,  and  indeed  even  their  very  presence.  Here  we  present  results  on  well- 
understood  molecular  assemblies,  which  exhibit  an  understood  classical  transport  behavior,  and 
which  can  be  used  as  a  control  for  eliminating  (or  understanding)  fabrication  variables.  Utilizing 
tunneling  spectroscopic  methods,  we  present  the  first  unambiguous  evidence  of  the  presence  of 
molecules  in  the  junction,  and  further  confirm  the  charge  transport  mechanism  obtained  by 
standard  current-voltage  characterizations. 

A  molecular  system  whose  structure  and  configuration  are  sufficiently  well-characterized 
such  that  it  can  serve  as  a  standard  is  the  extensively  studied  alkanethiol  (CH3(CH2)n-iSH)  self- 
assembled  monolayer  (SAM).9  This  system  is  useful  as  a  control  since  properly  prepared  SAMs 
form  single  van  der  Waals  crystals,9,10  and  presents  a  simple  classical  metal-insulator-metal  (M- 
I-M)  tunnel  junction  when  fabricated  between  metallic  contacts  due  to  the  large  HOMO-LUMO 
gap  (HOMO:  highest  occupied  molecular  orbital,  LUMO:  lowest  unoccupied  molecular  orbital) 
of  approximately  8  eV. 11-13 

Various  surface  analytical  tools  have  been  utilized  to  investigate  the  surface  and  bulk 
properties  of  the  alkanethiol  SAMs,  such  as  X-ray  photoelectron  spectroscopy,14  Fourier 
transform  infrared  spectroscopy  (FTIR),15  Raman  spectroscopy,16  scanning  tunneling  microscopy 
(STM)10  etc.  For  example,  studies  have  shown  that  the  bonding  of  the  thiolate  group  to  the  gold 
surface  is  strong  with  a  bonding  energy  of  -  1.7  eV.9  STM  topography  examinations  revealed 
that  alkanethiols  adopt  the  commensurate  crystalline  lattice  characterized  by  a  c(4  x  2) 

superlattice  of  a  (V 3  x  a/3)R30°.  10,17  FTIR  investigation  showed  that  the  orientation  of  the 

1 8 

alkanethiol  SAMs  on  Au(l  1 1)  surfaces  are  tilted  -  30°  from  the  surface  normal. 

Electronic  transport  through  alkanethiol  SAMs  have  also  been  characterized  by  STM,19'20 
conducting  atomic  force  microscopy,”  -  mercury-drop  junctions,  cross- wire  junctions,  and 

electrochemical  methods.30-32  These  investigations  are  exclusively  at  ambient  temperature  - 
clearly  useful  -  but  insufficient  for  an  unambiguous  claim  that  the  transport  mechanism  is 
tunneling  (of  course  expected,  assuming  that  the  Fermi  levels  of  the  contacts  lie  within  the  large 
HOMO-LUMO  gap).  However  in  the  absence  of  temperature-dependent  current-voltage  (I(V,T)) 
characteristics,  other  conduction  mechanisms  (such  as  thermionic,  hopping,  or  filamentary 
conduction)  cannot  be  excluded  complicate  the  analysis,  and  thus  such  a  claim  is  premature. 

Utilizing  a  nanometer  scale  device  structure  that  incorporates  alkanethiol  SAMs,  we 
demonstrate  devices  that  allow  I(V,T)  and  structure-dependent  measurements  '  with  results 
that  can  be  compared  with  accepted  theoretical  models  of  M-I-M  tunneling.  The  use  of  this 
fabrication  approach  is  not  special  in  any  way  (other  than  that  we  have  so  far  found  it  to  be 
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successful)  -  indeed  we  stress  that  any  successful  device  fabrication  method  should  yield  the 
results  described  below  if  one  is  characterizing  the  intrinsic  molecular  transport  properties. 

The  electronic  transport  is  further  investigated  with  the  technique  of  inelastic  electron 
tunneling  spectroscopy  (IETS).34  IETS  was  developed  in  the  1960’s  as  a  powerful  spectroscopic 
tool  to  study  the  vibrational  spectrum  of  organic  molecules  confined  inside  metal-oxide-metal 
junctions.  '  In  our  study  IETS  is  utilized  for  the  purpose  of  molecule  identification,  chemical 
bonding  and  conduction  mechanism  investigations  of  the  “control”  SAMs.  The  exclusive 
presence  of  well-known  vibrational  modes  of  the  alkanes  used  are  direct  evidence  of  the 
molecules  in  the  device  structure,  something  that  has  to  date  only  been  inferred  (with  good 
reason,  but  nonetheless  not  unambiguously).  The  vibrational  modes,  exclusively  identified  as 
alkanes  (as  well  as  contact  modes)  are  difficult  to  interpret  in  any  other  way  other  than  as 
components  in  the  active  region  of  the  device.  The  inelastic  tunneling  spectra  also  demonstrate 
that  electronic  tunneling  occurs  through  the  molecules,  confirming  the  conduction  mechanism 
obtained  by  I(V,T)  characterizations.  The  specific  spectral  lines  also  yield  intrinsic  linewidths 
that  may  give  insight  into  molecular  conformation,  and  may  prove  to  be  a  powerful  tool  in  future 
molecular  device  characterization. 

Electronic  transport  measurements  on  alkanethiol  SAMs  were  performed  using  a  device 
structure  similar  to  one  reported  previously.33  34'40-42  In  this  device,  as  illustrated  in  Figure  1(a) 
(not  drawn  to  scale  in  the  relative  thickness),  a  number  of  molecules  (~  several  thousands)  are 

sandwiched  between  two  metallic  contacts.  This 
technique  provides  a  stable  device  structure  and  makes 
cryogenic  measurements  possible.  The  device 
fabrication  starts  with  a  high  resistivity  silicon  wafer 
with  low  stress  Si3N4  film  deposited  on  both  sides  by 
low  pressure  chemical  vapor  deposition  (LPCVD).  By 
standard  photolithography  processing,  a  suspended 
Si3N4  membrane  (size  of  40  pm  x  40  pm  and  thickness 
of  ~  70  nm)  is  fabricated  on  the  topside  of  the  wafer. 
Subsequent  e-beam  lithography  and  reactive  ion  etching 
creates  a  single  pore  with  a  diameter  of  tens  of 
nanometers  through  the  membrane.  As  the  next  step, 

150  nm  gold  is  thermally  evaporated  onto  the  topside  of 
the  wafer  to  fill  the  pore  and  form  one  of  the  metallic 
contacts. 

Figure  1.  Schematics  of  a  nanometer  scale  device  used  in  this 
study  (not  drawn  to  scale  in  the  relative  thickness),  (a)  Top 
schematic  is  the  cross  section  of  a  silicon  wafer  with  a 
nanometer  scale  pore  etched  through  a  suspended  silicon 
nitride  membrane.  Middle  and  bottom  schematics  show  a 
Au/SAM/Au  junction  formed  in  the  pore  area,  (b)  The 
structures  of  octanethiol  and  octanedithiol  are  shown  as 
examples. 


(a) 


.SH 

HS 

octanedithiol 


The  device  is  then  transferred  into  a  molecular  solution  to  deposit  the  SAM  layer.  For  our 
experiments,  a  -  5  mM  alkanethiol  solution  is  prepared  by  adding  ~  1 0  pL  alkanethiols  into  10 
mL  ethanol.4'  The  deposition  is  done  in  solution  for  24  hours  inside  a  nitrogen  filled  glove  box 
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with  an  oxygen  level  of  less  than  100  ppm.  Three  alkanemonothiol  molecules  of  different 
molecular  lengths-octanethiol  [CH^CHVbSH;  denoted  as  C8,  for  the  number  of  alkyl  units], 
dodecanethiol  [CH3(CH2)nSH,  denoted  as  C 12],  and  hexadecanethiol  [CF^CH+hsSFI,  denoted 
as  Cl 6]  and  one  alkanedithiol  molecule-octanedithiol  [HSiCFfifiSH,  denoted  as  C8-dithiol] 
were  used  to  form  the  active  molecular  components.43  As  representative  examples,  the  chemical 
structures  of  octanethiol  and  octanedithiol  are  shown  in  Figure  1(b). 

In  order  to  statistically  determine  the  pore  size,  test  patterns  (arrays  of  pores)  were  created 
under  similar  fabrication  conditions.  Figure  2  shows  a  scanning  electron  microscope  (SEM) 
image  of  such  test  pattern  arrays.  This  indirect  measurement  of  device  size  is  done  since  SEM 
examination  of  the  actual  device  can  cause  hydrocarbon  contamination  of  the  device  and 
subsequent  contamination  of  the  monolayer.  From  regression  analysis  of  298  pores,  the  device 
sizes  of  the  C8,  C12,  C16,  and  C8-dithiol  samples  are  determined  as  50  ±  8,  45  ±  2,  45  ±  2,  and 
51+5  mn  in  diameters,  respectively.  A  more  ideal  (less  parasitic)  C8  sample  supercedes  that  of 
previous  reports,  and  derived  parameters  from  the  two  data  sets  agree  to  within  a  standard 

error.  We  will  use  these  device 
areas  as  the  effective  contact 
areas.  Although  one  could 
postulate  that  the  actual  area  of 
metal  that  contacts  the 
molecules  may  be  different, 
there  is  little  reason  to  propose  it 
would  be  different  as  a  function 
of  length  over  the  range  of 
alkanethiols  used,  and  at  most 
would  be  a  constant  systematic 
error. 

Figure  2.  A  scanning  electron 
microscope  image  of  a 
representative  array  of  pores  used 
to  calibrate  device  size.  The  scale 
bar  is  500  nm. 


The  sample  is  then  transferred  in  ambient  conditions  to  an  evaporator  that  has  a  cooling 
stage  to  deposit  the  opposing  Au  contact.  During  the  thermal  evaporation  (under  the  pressure  of 

o 

~  10"  Torr),  liquid  nitrogen  is  kept  flowing  through  the  cooling  stage  in  order  to  avoid  thermal 
damage  to  the  molecular  layer.33'44  This  technique  reduces  the  kinetic  energy  of  evaporated  Au 
atoms  at  the  surface  of  the  monolayer,  thus  preventing  Au  atoms  from  punching  through  the 
monolayer.  For  the  same  reason  the  evaporation  rate  is  kept  very  low.  For  the  first  10  nm  gold 
evaporated,  the  rate  is  less  than  0. 1  A/s.  Then  the  rate  is  increased  slowly  to  0.5  A/s  for  the  rest 
of  the  evaporation  and  a  total  of  200  nm  gold  is  deposited  to  form  the  contact. 

The  device  is  subsequently  packaged  and  loaded  into  a  low  temperature  cryostat.  The 
sample  temperature  is  varied  from  300  to  4.2  K  by  flowing  cryogen  vapor  onto  the  sample  (and 
thermometer)  using  a  closed  loop  temperature  controller.  Two-terminal  dc  I(V)  measurements 
are  performed  using  a  semiconductor  parameter  analyzer.  Inelastic  electron  tunneling  spectra  are 
obtained  via  a  standard  lock-in  second  harmonic  measurement  technique.  '  A  synthesized 
function  generator  is  used  to  provide  both  the  modulation  and  the  lock-in  reference  signal.  The 
second  harmonic  signal  (proportional  to  d2I/dV2)  is  directly  measured  using  a  lock-in  amplifier, 
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which  is  checked  to  be  consistent  with  a  numerical  derivative  of  the  first  hannonic  signal 
(proportional  to  dl/dV).  Various  modulation  amplitudes  and  frequencies  are  utilized  to  obtain  the 
spectra.  The  ac  modulation  is  added  to  a  dc  bias  using  operational  amplifier-based  custom 
circuitry.45 

In  Table  1,  possible  conduction  mechanisms  are  listed  with  their  characteristic  current, 
temperature-  and  voltage-dependencies46  (We  do  not  discuss  filamentary  tunneling  mechanisms, 
which  are  easier  to  categorize47).  Based  on  whether  thennal  activation  is  involved,  the 
conduction  mechanisms  fall  into  two  distinct  categories:  (i)  thennionic  or  hopping  conduction 
which  has  temperature-dependent  I(V)  behavior  and  (ii)  direct  tunneling  or  Fowler-Nordheim 
tunneling  which  does  not  have  temperature-dependent  I(V)  behavior.  For  example,  thennionic 
and  hopping  conductions  have  been  observed  for  4-thioacetylbiphenyl  SAMs40  and  1,4- 
phenelyene  diisocyanide  SAMs.41b  On  the  other  hand,  the  conduction  mechanism  is  expected  to 
be  tunneling  when  the  Fermi  levels  of  contacts  lie  within  the  large  HOMO-LUMO  gap  for  short 
length  molecule,  as  for  the  case  of  alkanethiol  molecular  system.1 113  Previous  work  on 

AO 

Langmuir-Blodgett  alkane  monolayers  exhibited  a  significant  impurity-dominated  transport 
component,  complicating  the  analysis.  I(V)  measurements  on  self-assembled  alkanethiol 
monolayers  have  also  been  reported;19'29'40  however  all  of  these  measurements  were  performed  at 
fixed  temperature  (300  K)  which  is  insufficient  to  prove  tunneling  as  the  dominant  mechanism. 


TABLE  1.  Possible  conduction  mechanisms.  Adapted  from  Ref.  46. 


Conduction 

Mechanism 
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Dependence 
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Dependence 

Direct 

Tunneling* 
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Hopping 

Conduction 
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*  This  characteristic  of  direct  tunneling  is  valid  for  the  low  bias  regime  [see  Eq.  (3a)]. 


To  describe  the  transport  through  a  molecular  system  having  HOMO  and  LUMO  energy 
levels,  one  of  the  applicable  models  is  the  Franz  two-band  model.50'53  This  model  provides  a 
non-parabolic  energy-momentum  E(k)  dispersion  relationship  by  considering  the  contributions 
of  both  the  HOMO  and  LUMO  energy  levels:50 


k2 


2  m '' 


E(  1  + 
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where  k  is  the  imaginary  part  of  wave  vector  of  electrons,  m*  is  the  electron  effective  mass, 
h  (=  27th)  is  Planck’s  constant,  E  is  the  electron  energy,  and  Eg  is  the  HOMO-LUMO  energy  gap. 
From  this  non-parabolic  E(k)  relationship,  the  effective  mass  of  the  electron  tunneling  through 
the  SAM  can  be  deduced  by  knowing  the  barrier  height  of  the  metal-SAM-metal  junction. 

When  the  Fermi  level  of  the  metal  is  aligned  close  enough  to  one  energy  level  (either 
HOMO  or  LUMO),  the  effect  of  the  other  distant  energy  level  on  the  tunneling  transport  is 
negligible,  and  the  widely  used  Simmons  model  4  is  an  excellent  approximation.55  Simmons 
model  expressed  the  tunneling  current  density  through  a  barrier  in  the  tunneling  regime  of  V  < 

®B/e  as25-54 


J  = 


y  Ak  2  h  d 1  j 
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eV 
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where  m  is  the  electron  mass,  d  is  the  barrier  width,  <DB  is  the  barrier  height,  and  V  is  the 
applied  bias.  For  molecular  systems,  the  Simmons  model  has  been  modified  with  a  parameter 
a."  ’  a  is  a  unitless  adjustable  parameter  that  is  introduced  to  provide  either  a  way  of  applying 
the  tunneling  model  of  a  rectangular  barrier  to  tunneling  through  a  nonrectangular  barrier,"'  or  an 
adjustment  to  account  for  the  effective  mass  (m*)  of  the  tunneling  electrons  through  a 
rectangular  barrier,  ’  ’  ’  or  both,  a  =  1  corresponds  to  the  case  for  a  rectangular  barrier  and 
bare  electron  mass.  By  fitting  individual  I(V)  data  using  Eq.  (2),  ®B  and  a  values  can  be 
obtained. 

Eq.  (2)  can  be  approximated  in  two  limits:  low  bias  and  high  bias  as  compared  with  the 
barrier  height  ®B.  For  the  low  bias  range,  Eq.  (2)  can  be  approximated  as54 
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(2  mOB) 


1/2 


e2a 


h  d 
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^f-a(v,rd 
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To  determine  the  high  bias  limit,  we  compare  the  relative  magnitudes  of  the  first  and  second 
exponential  terms  in  Eq.  (2).  At  high  bias,  the  first  term  is  dominant  and  thus  the  current  density 
can  be  approximated  as 


J 


d~ 


eV 


exp 


2(2/« ) 


1/2  f 

—  a 


- 


eV 


V'2 


d 


(3b) 


The  tunneling  currents  in  both  bias  regimes  are  exponentially  dependent  on  the  barrier  width 

d.  In  the  low  bias  regime  the  tunneling  current  density  is  J  oc  —  exp (~/3Qd)  ,  where  (30  is  bias- 

d 

independent  decay  coefficient: 


2(2/77)' 

n 

1 


-a 


(®.) 


1/2 


(4a) 


while  in  the  high  bias  regime,  J  oc  —  exp(-/?,,t/) ,  where  pv  is  bias-dependent  decay 
coefficient: 


d2 
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2(2777  ) 
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eV 
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eV 
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(4b) 


At  high  bias  Pv  decreases  as  f)ias  increases,  which  results  from  barrier  lowering  effect  due  to 
the  applied  bias. 


In  order  to  detennine  the  conduction  mechanism  of  self-assembled  alkanethiol  molecular 
systems  I(V)  measurements  in  a  sufficiently  wide  temperature  range  (300  to  80  K)  and  resolution 
(10  K)  were  performed.  Figure  3(a)  shows  a  representative  I(V,T)  characteristic  of  dodecanethiol 
(C12)  measured  with  the  device  structure  as  shown  in  Figure  1(a).  Positive  bias  corresponds  to 
electrons  injected  from  the  physisorbed  Au  contact  [bottom  contact  in  Figure  1(a)]  into  the 
molecules.  By  using  the  contact  area  of  45  ±  2  nm  in  diameter  determined  from  SEM  study,  a 
current  density  of  1,500  ±  200  A/cm  at  1.0  Volt  is  determined.  No  significant  temperature 
dependence  of  the  characteristics  (from  V  =  0  to  1.0  Volt)  is  observed  over  the  range  from  300  to 
80  K.  An  Arrhenius  plot  (ln(I)  versus  1/T)  of  this  is  shown  in  Figure  3(b),  exhibiting  little 
temperature  dependence  in  the  slopes  of  ln(I)  versus  1/T  at  different  bias  and  thus  indicating  the 
absence  of  thennal  activation.  Therefore,  we  conclude  that  the  conduction  mechanism  through 
alkanethiol  is  tunneling  contingent  on  demonstrating  a  correct  molecular  length  dependence.  The 
tunneling  through  alkanethiol  SAMs  has  been  assumed  as  “through-bond”  tunneling,  i.e.,  along 
the  tilted  molecular  chains  between  the  metal  contacts.”  ’  ’  ’  Based  on  the  applied  bias  as 
compared  with  the  barrier  height  (Ob),  the  tunneling  through  a  SAM  layer  can  be  categorized 
into  either  direct  (V  <  <J>B/e)  or  Fowler-Nordheim  (V  >  0B/e)  tunneling.  These  two  tunneling 
mechanisms  can  be  distinguished  due  to  their  distinct  voltage  dependencies  (see  Table  1). 
Analysis  of  ln(I/V2)  versus  1/V  [in  Figure  3(c)]  shows  no  significant  voltage  dependence, 
indicating  no  obvious  Fowler-Nordheim  transport  behavior  in  this  bias  range  (0  to  1.0  Volt)  and 
thus  detennining  that  the  barrier  height  is  larger  than  the  applied  bias,  i.e.,  Ob  >  1.0  eV.  This 
study  is  restricted  to  applied  biases  <1.0  Volt  and  the  transition  from  direct  to  Fowler-Nordheim 
tunneling  requires  higher  bias. 
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Figure  3.  (a)  Temperature- 
dependent  I(V) 
characteristics  of 
dodecanethiol  (C12).  I(V) 
data  at  temperatures  from 
300  to  80  K  with  20  K  steps 
are  plotted  on  a  log  scale,  (b) 
Arrhenius  plot  generated 
from  the  I(V)  data  in  (a),  at 
voltages  from  0.1  to  1.0  Volt 
with  0.1  Volt  steps,  (c)  Plot 
of  ln(I/V2)  versus  1/V  at 
selected  temperatures. 

The  importance  of 
variable  temperature 
measurements  to  validate 
tunneling  is  demonstrated 
in  Figure  4.  Here  the  I(V) 
of  an  octanethiol  (C8) 
device  is  shown  (Figure 
4(a)),  whose  I(V)  shape 
looks  very  similar  to 
Figure  3  (i.e.,  direct 
tunneling),  and  indeed  can 
be  fit  to  a  Simmons  model. 
However,  further  I(V,T) 
measurements  display  an 
obvious  temperature 
dependence  (Figure  4(b)), 
which  can  be  lit  well  to  a 
hopping  conduction  model 
(Table  1)  with  a  well- 
defined  activation  energy 
of  190  tneV  (Figure  4(c)). 
This  and  other  similar 
impurity-mediated 
transport  phenomena  (such 
as  Coulomb  blockade)  are 
observed  in  a  subset  of 
devices  and  is  indicative  of 
the  unintentional 
incorporation  of  a  trap  or 
defect  level  in  those 
devices.  This  study 
instead  focuses  on  devices 


which  do  not  show  any 

defect-mediated  transport  and  probes  the  intrinsic  behavior  of  the  molecular  layer. 
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Figure  4.  (a)  I(V)  characteristics 
of  an  octanethiol  (C8)  device  at 
270K.  (b)  Temperature 
dependence  of  the  device  from 
270K  to  18 OK  (in  10K 
increments),  (c)  Plot  of  ln(I/V) 
versus  1/T  at  various  voltages. 
The  activated  behavior  is  bias 
voltage  independent,  thus  the 
behavior  is  hopping  (in  this 
device)  due  to  incoiporation  of  a 
defect  of  energy  1 90  meV.  This 
class  of  device  is  not  suitable  for 
investigation  of  the  intrinsic 
transport  mechanism  in  the  SAM 
as  it  is  dominated  by  a  defect. 

Having  established 
tunneling  as  the  conduction 
mechanism  in  a  device,  we 
will  now  obtain  the  barrier 
height  by  comparing 
experimental  I(V)  data  with 
theoretical  calculations  from 
tunneling  models. 

From  the  modified 
Simmons  model  [Eq.  (2)]  by 
adjusting  two  parameters  Ob 
and  a,  a  nonlinear  least 
squares  fitting  can  be 
performed  to  fit  the  measured 
C12  I(V)  data  (calculation 
assuming  a  =  1  has  been 
previously  shown  not  to  fit 
I(V)  data  well  for  some 
alkanethiol  measurements  at 
fixed  temperature  (300  K)).25 
By  using  a  device  size  of  45 
nm  in  diameter,  the  best  fitting 
parameters  (minimizing  %  )  for 
the  room  temperature  C12  I(V) 
data  were  found  to  be  Ob  = 
1.42  ±  0.04  eV  and  a  =  0.65  ± 
0.01,  where  the  error  ranges  of 
Ob  and  a  are  dominated  by 
potential  device  size 
fluctuations  of  2  nm.  Likewise, 
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data  sets  were  obtained  and  fittings  were  done  for  octanethiol  (C8)  and  hexadecanethiol  (Cl 6), 
which  yielded  values  {Ob  =  1.83  ±  0.10  eV  and  a  =  0.61  ±0.01}  and  {®b  =  1.40  ±  0.03  eV,  a  = 
0.68  ±0.01},  respectively. 

Using  d>B  =  1.42  eV  and  a  =  0.65,  a  calculated  I(V)  for  C12  is  plotted  as  a  solid  curve  in  a 
linear  scale  [Figure  5(a)]  and  in  a  semi-log  scale  [Figure  5(b)].  A  calculated  I(V)  for  a  =  1  and 
Ob  =  0.65  eV  (which  gives  the  best  fit  at  low  bias  range)  is  shown  as  the  dashed  curve  in  the 
same  figure,  illustrating  that  with  a  =  1  only  limited  regions  of  the  I(V)  can  be  fit  (specifically 
here,  for  V  <  0.3  Volt).  For  the  case  of  a  rectangular  barrier,  the  a  parameter  fit  presented  above 
corresponds  to  an  effective  mass  m*  (  =  a  m)  of  0.42  m. 


Figure  5.  Measured  C12  I(V)  data  (circular  symbols)  is  compared  with  calculation  (solid  curve)  using  the 
optimum  fitting  parameters  of  ®B  =  1.42  eV  and  a  =  0.65.  The  calculated  I(V)  from  a  simple  rectangular 
model  (a  =  1 )  with  Ob  =  0.65  eV  is  also  shown  as  the  dashed  curve.  Current  is  plotted  (a)  on  linear  scale 
and  (b)  on  log  scale. 


In  order  to  investigate  the  dependency  of  the  Simmons  model  fitting  on  ®b  and  a,  a  fitting 
minimization  analysis  was  undertaken  on  the  individual  ®b  and  a  values  as  well  as  their  product 
form  of  a®B  “  in  Eq.  (4a).  A(®B,  a)  =  (  X  |  Iexp,v  -  Icai,v  |  )  was  calculated  and  plotted  where 
Iexp.v  is  the  experimental  current-voltage  values  and  Icai,v  is  calculated  using  Eq.  (2).  7,500 
different  {®B,  a}  pairs  were  used  in  the  fittings  with  ®B  ranging  from  1.0  to  2.5  eV  (0.01  eV 
increment)  and  a  from  0.5  to  F0  (0.01  increment).  Figure  6(a)  is  a  representative  contour  plot  of 
A(Ob,  a)  versus  ®b  and  a  values  generated  for  the  C 12  I(V)  data  where  darker  regions 
correspond  to  smaller  A(®B,  a)  and  various  shades  represent  half  order  of  magnitude  A(®B,  a) 
steps.  The  darker  regions  represent  better  fits  of  Eq.  (2)  to  the  measured  I(V)  data.  In  the  inset  in 
Figure  6(a)  one  can  see  there  is  a  range  of  possible  ®b  and  a  values  yielding  minimum  fitting 
parameters.  Although  the  tunneling  parameters  determined  from  the  previous  Simmons  tunneling 
fitting  {®b  =  1-42  eV  and  a  =  0.65}  lie  within  this  minimum  region  in  this  figure,  there  is  a 
distribution  of  other  possible  values. 
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Figure  6.  (a)  Contour  plot  of  A(Ob,  a) 
values  for  C12  nanopore  device  as  a 
function  of  ®B  and  a,  where  the  darker 
region  corresponds  to  a  better  fitting. 
Inset  shows  detailed  minimization 
fitting  regions,  (b)  A  plot  of  A(Ob,  a)  as 
a  function  of  a  0B12. 

A  plot  of  A(<I>b,  a)  versus 
a©B  ‘  for  the  same  device  reveals  a 
more  pronounced  dependence,  and  is 
shown  in  Figure  6(b).  This  plot 
indicates  the  fitting  to  the  Simmons 
model  sharply  depends  on  the 
product  of  a®B  •  For  this  plot  the 
A(®b,  ot)  is  minimized  at  a®B  “  of 
0.77  (eV)  corresponding  to  a  p0 
value  of  0.79  A'1  from  Eq.  (4a).  The 
C8  and  C 16  devices  showed  similar 
results,  indicating  the  Simmons 
tunneling  model  has  a  strong  a®B*  2 
dependence.  For  the  C8  device, 
although  ®B  obtained  from  the  fitting  is  a  little 
larger,  combined  a  and  ®B  gives  a  similar  p0 
value  within  the  error  range  as  the  C12  and  C16 
devices  (Table  2). 


TABLE  2.  Summary  of  alkanethiol  tunneling  parameters  in  this  study 


Molecules 

J  at  1  V  (A/cm2) 

®b  (eV) 

a 

m*  (m) 

Po(A4)* 

C8 

02 

31,000  ±  10,000 

1,500  ±200 

1.83  ±0.10 

1.42  ±0.04 

0.61  ±0.01 

0.65  ±0.01 

0.37 

0.42 

0.85  ±0.04 

0.79  ±  0.02 

06 

23  ±2 

1.40  ±0.03 

0.68  ±0.01 

0.46 

0.82  ±  0.02 

C8-dithiol 

93,000  ±  18,000 

1.20  ±0.03 

0.59  ±0.01 

0.35 

0.66  ±  0.02 
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*  p0  values  were  calculated  using  Eq.  (4a). 

Three  alkanethiols  of  different  molecular  length,  C8,  C12,  and  C16  were  investigated  to 
study  length-dependent  tunneling  behavior.  Figure  7  is  a  semi-log  plot  of  tunneling  current 
densities  multiplied  by  molecular  length  (Jd  at  low  bias  and  Jd2  at  high  bias)  as  a  function  of  the 
molecular  length  for  these  alkanethiols.  The  molecular  lengths  used  in  this  plot  are  13.3,  18.2, 
and  23.2  A  for  C8,  C12,  and  C16,  respectively.  Each  molecular  length  was  determined  by  adding 
an  Au-thiol  bonding  length  to  the  length  of  molecule  “.  Note  that  these  lengths  assume  through- 
bond  tunneling.  ’  ’  ’  The  high  and  low  bias  regimes  are  defined  somewhat  arbitrarily  by 
comparing  the  relative  magnitudes  of  the  first  and  second  exponential  terms  in  Eq.  (2).  Using  <bis 
=  1.42  eV  and  a  =  0.65  obtained  from  nonlinear  least  squares  fitting  of  the  C12  I(V)  data,  the 
second  term  becomes  less  than  ~  10  %  of  the  first  tenn  at  ~  0.5  Volt  that  is  chosen  as  the 
boundary  of  low  and  high  bias  ranges. 


Figure  7.  Log  plot  of  tunneling  current  densities  multiplied  by  molecular  length  d  at  low  bias  and  by  d2  at 
high  bias  (symbols)  versus  molecular  lengths.  The  lines  through  the  data  points  are  linear  fittings. 


As  seen  in  Figure  7,  the  tunneling  current  shows  exponential  dependence  on  molecular 
length,  which  is  consistent  with  the  Simmons  tunneling  model  [Eq.  (3)].  The  P  values  can  be 
determined  from  the  slope  at  each  bias  and  are  plotted  in  Figure  8.  The  error  bar  of  an  individual 
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P  value  in  this  plot  was  obtained  by  considering  both  the  device  size  uncertainties  and  the  linear 

fitting  errors. 


Figure  8.  Plot  of  p  versus 
bias  in  the  low  bias  range 
(square  symbols)  and  high 
bias  ranges  (circular 
symbols).  The  inset  shows  a 
plot  of  p2v  versus  bias  with 
a  linear  fitting. 

The  detennined  P 
values  are  almost 
independent  of  bias  in  the 
low  bias  range  (V  <  ~  0.5 
”  0.0  0.2  0.4  0.6  0.8  1.0  V),  and  an  average  P  of 

Y  (V)  0.77  ±  0.06  A'1  in  this 

region  (from  0  to  0.5  V) 
can  be  calculated  from 

Figure  8.  Table  3  is  a  summary  of  previously  reported  alkanethiol  transport  parameters  obtained 
by  different  techniques.  The  current  densities  (J)  listed  in  Table  3  are  for  C 12  mono  thiol  or 
dithiol  devices  at  1  V,  which  are  extrapolated  from  published  results  of  other  length  alkane 
molecules.  The  large  variation  of  J  of  these  reports  can  be  attributed  to  the  uncertainties  in  device 
contact  geometry  and  junction  area,  as  well  as  complicating  inelastic  or  defect  contributions. 

The  P  value  (0.77  ±  0.06  A'1  «  0.96  ±  0.08  per  methylene)  for  alkanethiols  reported  here  is 
comparable  to  previously  reported  values  as  summarized  in  Table  3.  This  P  value  agrees  with  the 
value  of  0.79  A'1  (p0)  calculated  via  Eq.  (4a)  from  fitting  individual  I(V)  characteristic  of  the  C12 
device.  The  calculated  Po  of  C8  and  C16  devices  also  have  similar  values,  as  summarized  in 
Table  2. 


TABLE  3.  Summary  of  alkanethiol  tunneling  characteristic  parameters 


Junction 

P  (A1) 

J  (A/cnr)  at  1  V 

Ob  (eV) 

Technique 

Ref. 

(bilayer)  monothiol 

0.87±0. 1 

25-200a) 

2.1e> 

Hg-j  unction 

25 

(bilayer)  monothiol 

0.71±0.08 

0.7-3.5a) 

Hg-junction 

27 

monothiol 

0.7910.01 

1500l200b) 

1.4e> 

Solid  M-I-M 

33 

monothiol 

1.2 

STM 

19 

dithiol 

0.810.08 

3.7-5  x  105c) 

512° 

STM 

20 

monothiol 

0.73-0.95 

1 100-1900d) 

2.2e> 

CAFM 

21 

monothiol 

0.64-0.8 

10-50d) 

2.3e> 

CAFM 

23 

dithiol 

0.4610.02 

3-6  x  105c) 

1.3-1.5e) 

CAFM 

24 
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monothiol 

1.37±0.03 

1.8° 

Tuning  fork  AFM 

49 

monothiol 

0.97±0.04 

Electrochemical 

30 

monothiol 

0.85 

Electrochemical 

31 

monothiol 

0.91±0.08 

Electrochemical 

32 

monothiol 

0.76  2xl04(at  0.1  V)c) 

1.3-3.48’ 

Theory 

58 

monothiol 

0.76 

Theory 

59 

monothiol 

0.79 

Theory 

56 

Note: 

Some  decay  coefficients  P  were  converted  into  the  unit  of  A"1 2  from  the  unit  of  per  methylene. 

The  junction  areas  were  estimated  by  optical  microscope3',  SEMb),  assuming  single  molecule3’,  and  Flertzian  contact 
theoryd). 

Current  densities  (J)  for  C12  monothiol  or  dithiol  at  1  V  are  extrapolated  from  published  results  for  other  length 
molecules  by  using  conductance  oc  exp(-(l  d)  relationship. 

Barrier  height  <DB  values  were  obtained  from  Simmons  equation3’,  bias-dependence  of  (V1.  and  a  theoretical 
calculation8’. 

2 

According  to  Eq.  (4b),  P  y  depends  on  bias  V  linearly  in  the  high  bias  range.  The  inset  in 
Figure  8  is  a  plot  of  P'y  versus  V  in  this  range  (0.5  to  1.0  Volt)  along  with  linear  fitting  of  the 
data.  From  this  fitting,  <bB  =  1.35  ±  0.20  eV  and  a  =  0.66  ±  0.04  were  obtained  from  the  intercept 
and  the  slope,  respectively,  consistent  with  the  values  {<f>B  =  1.42  eV  and  a  =  0.65}  obtained 
from  the  nonlinear  least  squares  fitting  in  the  previous  section. 

P  values  for  alkanethiols  obtained  by  various  experimental  techniques  have  previously  been 
reported  and  are  summarized  in  Table  3.19'33’49  In  order  to  compare  with  these  reported  P  values, 
we  also  perfonned  length-dependent  analysis  on  our  experimental  data  according  to  the  generally 

.  ■  •  20-28,33 

used  equation: 


G  =  G0  exp(-/?  d )  (5) 

This  gives  a  P  value  from  0.84  to  0.73  A'1  in  the  bias  range  from  0.1  to  1.0  volt,  which  is 
comparable  to  results  reported  previously.  For  example,  Holmlin,  et.  al,  reported  a  P  value  of 
0.87  A'1  by  mercury  drop  experiments,25  and  Wold,  et.  al,  have  reported  P  of  0.94  A'1  and  Cui, 
et.  al,  reported  P  of  0.64  A'1  for  various  alkanethiols  by  using  a  conducting  atomic  force 
microscope  technique.'  ’  These  reported  P  were  treated  as  bias-independent  quantities,  contrary 
to  the  results  reported  here  and  that  observed  in  a  slightly  different  alkane  system  (ligand- 
encapsulated  nanoparticle/alkane-dithiol  molecules).24  We  also  caution  against  the  use  of 
parameters  that  have  not  been  checked  with  a  temperature-dependent  analysis,  since  small  non¬ 
tunneling  components  can  dramatically  affect  derived  values  of  p. 

We  have  analyzed  our  experimental  data  using  a  Franz  two-band  model.50'53  Since  there  is 
no  reliable  experimental  data  on  the  Fermi  level  alignment  in  these  metal-SAM-metal  systems, 
Ob  and  m*  are  treated  as  adjustable  parameters.  We  performed  a  least  squares  fit  on  our  data 
with  the  Franz  non-parabolic  E(k)  relationship  [Eq.  (1)]  using  an  alkanethiol  HOMO-LUMO  gap 
of  8  eV.  ’  Figure  9  shows  the  resultant  E(k)  relationship  and  the  corresponding  energy  band 
diagrams.  The  zero  of  energy  in  this  plot  was  chosen  as  the  FUMO  energy.  The  best  fitting 


15 


Final  report  FA9550-05-1-0395 

2 

parameters  obtained  by  minimizing  %  were  Ob  =  1.49  ±  0.5 1  eV  and  m*  =  0.43  ±  0. 15  m,  where 
the  error  ranges  of  Ob  and  m*  are  dominated  by  the  error  fluctuations  off)  [k  =  -  (p/2)”].  Both 
electron  tunneling  near  the  LUMO  and  hole  tunneling  near  the  HOMO  can  be  described  by  these 

parameters.  Ob  =  1 .49  eV 
indicates  that  the  Fermi  level  is 
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aligned  close  to  one  energy  level 
in  either  case,  therefore  the 
Simmons  model  is  a  valid 
approximation.  The  Ob  and  m* 
values  obtained  here  are  in 
reasonable  agreement  with  the 
previous  results  obtained  from  the 
Simmons  model. 

Figure  9.  E(k)  relationship  (symbols) 
generated  from  the  length-dependent 
measurement  data  for  alkanethiols. 
Solid  and  open  symbols  correspond 
to  electron  and  hole  tunneling, 
respectively.  The  insets  show  the 
corresponding  energy  band  diagrams. 
The  solid  curve  is  the  Franz  two- 
band  expression  for  m*  =  0.43  m. 


Electronic  transport  through  alkanethiol  SAMs  is  further  investigated  with  the  technique  of 
inelastic  electron  tunneling  spectroscopy,34  such  as  the  works  of  1966  by  Jaklevic  and  Lambe 
who  studied  the  conductance  of  a  tunnel  junctions  with  encased  organic  molecules.  Since  then 
it  has  become  a  powerful  spectroscopic  tool  for  chemical  identification,  chemical  bonding 
investigation,  and  surface  chemistry  and  physics  studies.  In  an  inelastic  tunneling  process  the 
electron  loses  energy  to  a  localized  vibrational  mode  with  a  frequency  v  when  the  applied  bias 
satisfies  the  condition  of  eV  =  hv.  As  a  result,  an  additional  tunneling  channel  is  opened  for  the 
electron,  resulting  in  an  increase  in  the  total  current  at  the  applied  bias  corresponding  to  the 
vibrational  mode  energy.  Typically  only  a  small  fraction  of  tunneling  electrons  are  involved  in 
the  inelastic  tunneling  process  (detennined  by  the  electron  -  vibronic  mode  coupling 
coefficient),  resulting  in  a  small  conductance  change,  which  is  commonly  measured  in  the 
second  harmonics  of  a  phase-sensitive  detector  that  yields  the  characteristic  frequencies  of  the 

36  38 

corresponding  vibrational  modes  as  well  as  other  information/ 
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I(V,T)  measurements  and 
additional  IETS  studies  have  been 
performed  on  an  octanedithiol  (C8- 
dithiol)  SAM  using  the  aforementioned 
device  structure  shown  in  Figure  1(a).34 
Figure  10(a)  is  the  I(V,T)  data  for  this 
device  obtained  from  300  to  4.2  K.  An 
Arrhenius  plot  shown  in  Figure  10(b) 
exhibits  little  temperature  dependence, 
verifying  that  tunneling  is  the  main 
transport  mechanism  for  C8-dithiol 
SAM.  This  result  is  in  good  agreement 
with  the  tunneling  transport 
characteristics  observed  previously. 
Figure  10(c)  shows  the  room 
temperature  I(V)  measurement  result. 
Using  a  junction  area  of  5 1  ±  5  nm  in 
diameter  (obtained  from  statistical 
studies  of  the  nanopore  size  with  SEM), 
a  current  density  of  (9.3  ±  1.8)  x  104 
A/cm2  at  1.0  Volt  is  calculated.  As  a 
comparison,  the  current  density  of  (3.1 
±  1.0)  x  104A/cm2  at  1.0  Volt  was 
observed  for  C8  monothiol  SAM. 

Using  the  modified  Simmons  model 
[Eq.  (2)],  the  transport  parameters  of 
©b  =  1.20  ±  0.03  eV  and  a  =  0.59  ± 
0.01  (m*  =  0.34  m)  were  obtained  for 
this  C8-dithiol  SAM. 


Figure  10.  (a)  I(V,T)  characteristics  of  C8 
dithiol  SAM  at  selected  temperatures  (4.2, 
50,  100,  150,  200,  250,  and  290  K).  (b) 
Arrhenius  plot  generated  from  the  data  in 
(a),  at  voltages  from  0.1  to  0.5  Volt  with 
0.05  Volt  steps,  (c)  Measured  C8-dithiol 
I(V)  data  at  room  temperature  (circular 
symbols)  is  compared  with  calculation 
(solid  curve)  using  the  optimum  fitting 
parameters  of  Ob  =  F20  eV  and  a  =  0.59. 


Figure  1 1  shows  the  IETS  spectrum  of  the  same  C8-dithiol  SAM  device  obtained  at  T  =  4.2 
K.  An  ac  modulation  of  8.7  mV  (rms  value)  at  a  frequency  of  503  Hz  was  applied  to  the  sample 
to  acquire  the  second  harmonic  signals.  The  spectra  are  stable  and  repeatable  upon  successive 
bias  sweeps.  The  spectrum  at  4.2  K  is  characterized  by  three  pronounced  peaks  in  the  0  to  200 
mV  region  at  33,  133,  and  158  mV.  From  comparison  with  previously  reported  infrared  (IR), 
Raman,  and  high  resolution  electron  energy  loss  (HREEF)  spectra  of  SAM  covered  gold  surfaces 
(Table  4),  these  three  peaks  are  assigned  to  v(Au-S),  v(C-C),  and  yw(CH2)  modes  of  a  surface 
bound  alkanethiolate. 60-63  The  absence  of  a  strong  v(S-H)  signal  at  -329  mV  suggests  that  most 
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of  the  thiol  groups  have  reacted  with  the  gold  bottom  and  top  contacts.  Peaks  are  also 
reproducibly  observed  at  80,  107  mV,  and  186  mV.  They  correspond  to  v(C-S),  5r(CH2),  and 
8s(CH2)  modes.  The  stretching  mode  of  the  CH2  groups,  v(CH2),  appears  as  a  shoulder  at  357 
rneV.  The  peak  at  15  mV  is  due  to  vibrations  from  either  Si,  Au,  or  §(C-C-C).64  We  note  that  all 
alkanethiolate  peaks  without  exception  or  omission  occur  in  the  spectra.  Peaks  at  58,  257,  277, 
and  302,  as  well  as  above  375  mV  are  likely  to  originate  from  Si-H  and  N-H  vibrations  related  to 
the  silicon  nitride  membrane,643'65  which  forms  the  SAM  encasement.  To  the  best  of  our 
knowledge  alkanethiols  have  no  vibrational  signatures  in  these  regions.  Measurement  of  the 
background  spectrum  of  an  “empty”  nanopore  device  with  only  gold  contacts  to  obtain 
background  contributions  from  ShN4  is  hampered  by  either  too  low  (open  circuit)  or  too  high 
(short  circuit)  currents  in  such  a  device.  Similar  IETS  result  has  also  been  obtained  using  a 
different  test  structure  recently.66 


Figure  11.  Inelastic  electron  tunneling  spectrum  of  C8  dithiol  SAM  obtained  from  lock-in  second 
harmonic  measurement  with  an  ac  modulation  of  8.7  mV  (rms  value)  at  a  frequency  of  503  Flz  (T  =  4.2 
K).  Peaks  labeled  with  Si  are  most  probably  background  due  to  the  encasing  ShN4. 


TABLE  4:  Summary*  of  the  major  vibrational  modes  of  alkanethiolates.  Taken  from  Ref.  60-62. 


Modes 

Methods  Wavenumber  (cm4)  (meV) 

v(Au-S)  HREELS61  225  28 

v(C-S)  Raman60  641  79 
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180 
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360 
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*  There  is  a  vast  amount  of  literature  with  spectroscopic  assignments  for  alkanethiols.  The 
references  given  are  representative  for  IR,  Raman,  and  HREELS  assignments. 


Although  there  are  no  selection  rules  in  IETS  as  there  are  in  IR  and  Raman  spectroscopy, 
certain  selection  preferences  have  been  established.  According  to  the  IETS  theory,  molecular 
vibrations  with  net  dipole  moments  perpendicular  to  the  interface  of  the  tunneling  junction  have 
stronger  peak  intensities  than  vibrations  with  net  dipole  moments  parallel  to  the  interface  (for 
dipoles  close  to  the  electrodes).  Thus  vibrations  perpendicular  to  the  electrode  interface,  i.e., 
v(Au-S),  v(C-S),  v(C-C),  and  yw(CH2),  dominate  the  IETS  spectrum  while  modes  parallel  to  the 
interface,  i.e.,  5r>s(CH2)  and  v(CH2),  are  weak,  as  clearly  shown  in  Figure  11. 

In  order  to  verify  that  the  observed  spectra  are  indeed  valid  IETS  data,  the  peak  width 
broadening  was  examined  as  a  function  of  temperature  and  modulation  voltage.  IETS  was 
performed  with  different  ac  modulations  at  a  fixed  temperature,  and  at  different  temperatures 
with  a  fixed  ac  modulation.  Figure  12(a)  shows  the  modulation  dependence  of  the  IETS  spectra 
obtained  at  4.2  K,  and  Figure  12(b)  shows  the  modulation  broadening  of  the  C-C  stretching 
mode  at  133  meV.  The  circular  symbols  are  the  full  widths  at  half  maximum  (FWHMs)  of  the 
experimental  peak  at  T  =  4.2  K  with  various  modulation  voltages.  A  Gaussian  distribution 
function  was  utilized  to  obtain  a  FWHM  and  the  error  range.  The  square  symbols  are  calculated 
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FWHM  values  (Wthe0reticai)  taking  into  account  both  a  finite  temperature  effect  (Wthermai  ~  5.4 
kBT)36  and  a  finite  voltage  modulation  effect  (Wmoduiation  ~  1.7  Vac  rms)-69  These  two  broadening 
contributions  add  as  the  squares:  W2theoreticai  =  W2thermai  +  W2moduiation.  The  agreement  is  excellent 
over  most  of  the  modulation  range,  but  we  note  a  saturation  of  the  linewidth  at  low  modulation 
bias  indicating  the  influence  of  a  non-negligible  intrinsic  linewidth.  Taking  into  account  the 
known  thermal  and  modulation  broadenings,  and  including  the  intrinsic  linewidth  (Wi)70  as  a 
fitting  parameter,  the  measured  peak  width  (Wexp)  is  given  by 

Wexp  =  ylWI2  +  W thermal  +  W modulation  (6) 

Wi  can  be  determined  by  using  a  nonlinear  least  squares  fit  to  the  ac  modulation  data  (Figure 
12)  with  Eq.  (6),  giving  an  intrinsic  linewidth  of  3.73  ±  0.98  meV  for  this  line.  This  is  shown 
(with  the  error  range)  in  Figure  12(b)  as  a  shaded  bar,  including  the  thermal  contribution. 


0  1000  2000  3000  4000  cm 


11.6mV 

10.2mV 

8.7mV 

7.3mV 

5.8mV 

2.9mV 


T  =  4.2  K 


V(V) 


Figure  12.  (a)  Modulation  dependence  of  IETS 
spectra  obtained  at  4.2  K.  (b)  Line  (C-C 
stretching  mode)  broadening  as  a  function  of  ac 
modulation.  The  circular  symbols  are 
experimental  FWHMs  and  the  square  symbols 
are  theoretical  calculations  considering  both 

modulation  and  thermal  contributions.  The  shaded  bar  denotes  the  expected  saturation  due  to  the  derived 
intrinsic  linewidth  (including  a  5.4  kBT  thermal  contribution)  of  3.73  ±  0.98  meV. 
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We  can  independently  check  the  thermal  broadening  of  the  line  at  fixed  modulation.  Figure 
13(a)  shows  the  temperature  dependence  of  the  IETS  spectra  obtained  with  an  ac  modulation  of 
8.7  mV  (rms  value).  In  Figure  13(b)  the  circular  symbols  (and  corresponding  error  bars)  are 
experimental  FWFIM  values  of  the  C-C  stretching  mode  from  Figure  13(a),  determined  by  a 
Gaussian  fit  (and  error  of  the  fit)  to  the  experimental  lineshape.  For  simplicity  we  have  only 
considered  Gaussian  lineshapes68  resulting  in  increased  error  bars  for  the  lower  temperature 
range  due  to  an  asymmetric  lineshape.  The  square  symbols  are  theoretical  calculations 
considering  thermal  broadening,  modulation  broadening,  and  the  intrinsic  linewidth  determined 
above.  The  error  ranges  of  the  calculation  (due  to  the  intrinsic  linewidth  error)  are  approximately 
the  size  of  the  data  points.  The  agreement  between  theory  and  experiment  is  very  good,  spanning 
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a  temperature  range  from  below  (x  0.5)  to  above  (xlO)  the  thermally  broadened  intrinsic 
linewidth.  This  linewidth  should  be  a  sensitive  test  to  compare  to  theoretical  models  of 
transmission  probabilities.71 
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Figure  13.  (a)  Temperature  dependence  of  IETS 
spectra  obtained  at  a  fixed  ac  modulation  of  8.7 
mV  (rms  value),  (b)  Line  (C-C  stretching  mode) 
broadening  as  a  function  of  temperature.  The 
circular  symbols  are  experimental  FWHMs  and 
the  square  symbols  are  theoretical  calculations 
considering  thermal  broadening,  modulation 
broadening,  and  the  intrinsic  linewidth.  . 


Similar  intrinsic  linewidths  have  been  determined  for  the  Au-S  stretching  mode  (33  meV) 
and  the  CH2  wagging  mode  (158  meV),  and  the  results  are  shown  in  Figure.  14.  For  the  Au-S 
stretching  mode,  the  deviation  of  experimental  data  from  calculated  values  (thermal  and 
modulation  width  only)  is  little  (Fig.  14(a)),  indicating  that  its  intrinsic  linewidth  is  small.  A 

linewidth  upper  limit  of  1 .69  meV  is  determined  for 
this  vibrational  mode.  For  the  CFI2  wagging  mode,  a 
nonlinear  least  squares  fit  to  Eq.  (6)  (solid  curve  in 
Fig.  14(b))  gave  an  intrinsic  linewidth  of  13.5  ±  2.4 
meV.  The  linewidths  and  their  variation  throughout 
the  molecule  are  potentially  due  to  inhomogeneous 
configuration  of  the  molecular  constituents,  and  a 
more  detailed  understanding  may  give  detailed 
structural  information  of  these  device  structures. 


Figure.  14.  Line  broadenings  as  function  of  AC  modulation 
for  IETS  spectra  obtained  at  4.2  K  for  (a)  the  Au-S 
stretching  mode  and  (b)  the  CH2  wagging  mode.  The 
circular  symbols  are  experimental  FWHMs  and  the  square 
symbols  are  theoretical  calculations  considering  modulation 
and  thermal  contributions.  Nonlinear  least  squares  fitting  to 
determine  intrinsic  linewidth  is  shown  as  the  solid  curve  in 
(b).  The  intrinsic  linewidths  obtained  for  the  Au-S 
stretching  mode  is  <  1.69  meV,  and  for  the  CH2  wagging 
mode  is  13.5  ±  2.4  meV. 
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